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ABSTRACT 

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. 
Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the 
Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. 
Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable 
consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are 
caused by magnetic field rearrangements. 

Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art 
microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected 
only in the crust of neutron stars, where we focus our study. 

Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different 
stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ~ 3) or low 
field neutron stars that accreted fallback matter during the supernova explosion (systematically n < 3); in a second stage, the evolution 
is governed by almost pure Ohmic field decay, and a braking index n > 3 is expected; in the third stage, at late times, when the interior 
temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of high absolute 
value and both positive and negative signs. 

Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the 
braking index, from temporal variations of the magnetic field. Models with strong (> 10 14 G) multipolar or toroidal components, even 
with a weak (~ 10 12 G) dipolar field are consistent with the observed trend of the timing properties. 

Key words. pulsars:general - stars: neutron - stars: magnetic fields - stars: evolution 



1. Introduction. 

The spin-down of an isolated neutron star (NS) on secular 
timescales is mainly caused by rotational energy losses due to 
electromagnetic radiation, pulsar winds, or gravitational radia- 
tion. A measurable quantity closely related to the rotational evo- 
lution of pulsars is the braking index n, defined by assuming 
that the star spins down according to a power law Cl = -KQ.", 
where £2 is the NS angular velocity. For each of these dominant 
rotational energy loss mechanisms, the proportionality constant 
K hides different dependences on the star radius, moment of in- 
ertia, magnetic field strength, and angle between rotation and 
magnetic axis. If all these quantities are constant in time, the 
magneto-dipole spindown mechanism predicts a braking index 
n — 3, but variations in time of any of these quantities may cause 
departures from this canonical value. Unfortunately, the accurate 
determination of the second derivative of the frequency, needed 
to estimate the braking index, is not always possible because it 
is affected by glitches and other short-term timing irregularities. 

At present, eight pulsars have sufficiently steady rota- 
tions that stable values of their brakin g index are gen- 
erally accepted (iLvneetalJ 19931 [19%!: Middleditch et al.l 
2006 : iLivingstone et al.1 120071 l201lt IWeltevrede et al] l201lt 
Espinoza et al.ll201ll) . among which the most recent case is PSR 
J1734-3333, which h as n = 0.9 + 0.2 significantly below 3 
(lEspinoza et alJ20Tlh . All these cases show the same trend: they 



are all young pulsars (Vela is the oldest among these pulsars with 
t » 12 kyr) and they all have n < 3. We note that gravitational 
wave emission predicts n — 5, but it is only efficient during the 
first minutes or hours of a NS life, when rotation is sufficiently 
fast and the mass q uadrupole moment large enough (see e.g., 
lHaskell et al.1 J2006) and references therein). It certainly does not 
contribute significantly to the timing noise of pulsars older than 
100 yr. 

Selecti ng a sample of 127 pu lsars from the ATNF Pulsar 
Catalogue dManchester et al.l|2005l) . for which the quoted errors 
in the se cond derivative of th e spin frequency (v) are less than ten 
per cent. lUrama et al.ld20 06) found a strong correlation of v with 
v, independent of the sign of v. They suggested that this trend 
can be accounted for by small stochastic deviations in the spin- 
down torque that are directly proportional (in magnitude) to the 
spin-down torque itself. Another point discussed in the literature 
is that some of the old pulsars (> 10 yr) have braking indices 
with absolute values exceeding \n\ - 10 4 . The occurrence of very 
high braking indices of both signs has been considered in the 
context of internal frictional instabilities occurring between the 
crust and the superfluid, almost i ndependently of the evolution o f 
the neutron star magnetic field (Shibazaki & Mochizuki 19951) . 
However, this applies only for old neutron stars (t > 2 x 10 7 
yr) and appears as extrem ely short term events oscilla ting about 
the canonical value n — 3. Barsukov & Tsvgan (2010) proposed 
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another explanation for the observed distribution with very high 
positive and negative braking indices by studying the effect of 
nondipolar magnetic field components and neutron star preces- 
sion on magnetospheric electric current losses. These large n 
should be observable over relatively long periods of 10 3 -10 4 yr. 
Another possibility that can explain the observed variability of 
braking ind ices is the time- evolution of conductivity in the mag- 
netosphere (ILi et al.ll2012bl) . which also has implications for the 
spin-down o f intermittent pulsars and subpulse drift phenomena 
(lLvnell2009l) . 

In this paper we focus on the imprint that the time-evolution 
of the internal magnetic field has on the timing noise. Our goal is 
to estimate the contribution of this variability (expected in realis- 
tic evolutionary models) to the braking index at different epochs. 
In particular, we discuss the possible reemergence of a magnetic 
field initially submerged by hypercritical accretion during the 
supernova explosion, and the evolution of the crustal magnetic 
field under the combined action of Ohmic dissipation and Hall 
drift during the first few million years of a NS life. 

2. Basic equations. 

The spin-down behavior of a rotating neutron star is governed 
by the energy balance equation relating the loss of rotational 
energy, E ml = 7Q 2 /2, where / is the moment of inertia of the 
neutron star and Q the angular velocity, to the energy loss rate 
E, given by magneto-dipole radiation, wind, gravitational radia- 
tion, or others: 



E mt = IQ.D. * E. 

The standard way to define the braking index is 
_ QQ _ vv _ PP 



(1) 



n 2 



p 2 



where v is the spin frequency, and we denote by P = 1/v the 
period. Under the usual assumption of a power-law rotational 
evolution 



Q = -KQ? . 



(3) 



the power-law index coincides with the braking index. 

For the particular case of a rotating dipole in vacuum, the 
well-known Larmor formula returns 



E = 



BlR 6 Q. A 



3c 3 



■ sin a 



(4) 



where Bo is the dipolar component of the surface magnetic field 
at the pole, R denotes the neutron star radius, a is the angle be- 
tween the rotational and the magnetic axis, and c is the speed of 
light. We neglect the contribution to the spindown of higher or- 
der multipoles, since they have significantly shorter "lever arms" 
than the dipolar mode. 

Combining Eqs. ([1) and ©, one obtains 



Q = -KCl 3 , 



(5) 



where K = faB^R 6 /31c 3 . For simplicity, but without loss of gen- 
erality, we have omitted the geometric factor that carries the 
dependence on the inclination angle, and we assume hereafter 
that it is a correction f a factorized in K, which carries infor- 
mation about the particular physical process that governs en- 
ergy losses. The functional dependence of K for all other mag- 
netic processes is the same oc B 2 £2 4 B 6 . Differences in the radi- 
ation mechanism are included in the factor f a . While magneto - 
spheric current losses scale as cos 2 a (Beskin & Nokhrina 2007), 



magneto-dipole radiation losses scale as sin a for vacuum o r 
|(1 + sin 2 a) for force-free magnetospheres (Spitkovskv 2006). 
The most recen t resistive solutions for pulsar magnetospheres 
(ILi et al.ll2012bl) fit the spin-down luminosity with a prefactor of 
the order of unity that also depends on the maximum potential 
drop along field lines in the corotating frame. 
Eq. <(5j can also be cast in the usual form 



PP = KB%(f) 



(6) 



where for a standard neutron star (R = 10 6 cm, / = 10 45 gcm 2 ) the 
constant 7C « 10~ 39 cm s 3 g _1 . Here we have explicitly written 
Bo(t) to show that the magnetic field is the only quantity that we 
allow to vary with time. 

Now we consider what happens when a secular variation of 
the magnetic field is allowed. Deriving P from Eq. (0, we have 



P = <KBP 2 (2B {) P - B P) , 



(7) 



and the braking index in a general case of a time-dependent mag- 
netic field can be simply expressed as follows 



n = 3 _A cS 3-4^, 
Bo T B 



(8) 



where t c = P/2P is the characteristic age, and we defined the 
magnetic field evolution timescale tb = Bo/ Bo. This last equa- 
tion shows that any variation of the dipolar surface magnetic 
field strength results in a deviation from the n — 3 standard 
value, which is obviously recovered for a constant magnetic field 
(Bo = 0). For an increasing Bo we will always obtain n < 3, 
while n > 3 is the signature of a decreasing Bq. 



(2) 3. Magnetic field evolution scenarios. 



The precise calculation of the time variation of Bo requires solv- 
ing the proble m of the coupled magneto-thermal e volution of 
a neutron star dPons et al . 2009: IVigan6 et al.l [2012). However, 
before showing results from simulations, one can qualitatively 
analyze of the most important physical process at different ages. 
We now discuss three possible scenarios. 

3.1. Amplification or reemergence of the dipolar surface 
magnetic field. 

The generic n < 3 observed without exception for young 
pulsars can be caused by the rediffusion of the magnetic 
field submerg ed in the crust during the supernova fall- 
back e pisode ([Young & Chanmugaml [l 995 ; Muslimov & Page 
1 1995b I Geppert et al.l Il999t) . The submergence of the mag- 
netic field has also been studied in recent MHD simulations, 
with interesting implications for gravitation al wave emission 
dVigelius & Melatosl 12009b IWette et all l2010l) . The rediffusion 
timescale depends essentially on the total amount of accreted 
matter. For total accreted masses between 10" 4 -10" 3 M o , screen- 
ing currents are dissipated on T rediff ~ 10 3 -10 4 yr, as recen t cal- 
culations have confirmed (lHdl20 1 lb l\^"gano & Ponsll2012l) . The 
same idea has been applied to the reemergence o f magnetic fields 
in acc reting binary systems with similar results (Cumming et al.l 
120011) . and evidence for fast magnetic field evolution in an ac- 
creting millisecond pulsar has been reported (Pat runoll2012l) . 

Alternative mechanisms are time variat ions in the angle be- 
tween the rotational a nd magnetic axes ([Link & Eps teinlll997b 
iRuderman etaTI 1 19981) . or the thermoelectric field generation 
that may proceed in the crust and envelope of young pulsars if 
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a sufficiently strong temperature gradient is present (lUrpin et al.1 
l!986UWiebicke & Gepperll 996). This process is limited by the 
condition that the surface te mperature of the neutron star should 
not be lower than 3 x 10 6 K dWiebicke & Gepperll996l) . which 
corresponds to about 1000 yr in the standard cooling scenario. 
In our simulations we included the reemergence of a screened 
magnetic field, but our present version of the code does not in- 
clude magnetic field generation by thermoelectric effect, or other 
possible mechanisms such as magnetic flux expulsion from the 
superconducting core. 

3.2. Ohmic decay. 

Ohmic dissipation of the magnetic field has been thoroughly 
investigated in the past by many grou ps (lHaensel et al.l l 1990; 
iPage et al.|[2QQ0i: iTauris & K onar 2001). The usual definition of 
the Ohmic decay time is 



modes are expected. The temporal evolution of the polar surface 
magnetic field can then be approximated by 



Ano-L 1 



1"Ohn 



(9) 



where L denotes the scale-length of magnetic field variations 
and cr the electric conductivity. In a neutron star crust, cr is 
dominated by electronic transport and depends on the electron 
density, the crustal temperature, and the impurity concentration 
within the crust. Since the electron density varies over about 
four orders of magnitude in the crust and the temperature de- 
creases by about two to three orders of magnitude in a pulsar's 
lifetime, the electric conductivity may vary b oth in space and 
time b y many orders of magnitude (see Fig. 1 in lPons & Geppertl 
d2007l) ). Therefore, to assume a uniform Ohmic decay time rohm, 
independent of the location of currents and the pulsar age is, gen- 
erally speaking, quite misleading: Ohmic decay in a neutron star 
crust cann o t be d escribed by a single exponential law (see e.g. 
IPage et ail j2000t) for a qualitative discussion). Since the con- 
ductivity increases with decreasing temperature, Ohmic diffu- 
sion becomes an increasingly slow process in older pulsars, be- 
yond the transition from neutrino dominated to photon cooling 
at > 10 5 yr. The only general fact that we can expect is that, for 
low-field middle-aged neutron stars (10 4 -10 5 yr), this is expected 
to be the dominant process, thus resulting in n > 3. 

3.3. Hall drift oscillatory modes. 

The Hall drift, through its nonlinear dependence on the magnetic 
field, has an influence on the magnetic field evolution either for 
magnetar conditions (B > 10 14 G), or alternatively when the 
electrical resistivity becomes very low. Under certain circum- 
stances the Hall drift can drain magnetic energy out of the dipo- 
lar mode and redistribute it into smaller scale ones. This some- 
times causes an apparently rapid decrease of Bo- But the oppo- 
site effect may also occur, i.e., pumping energy from an internal 
strong toroidal field to the dipolar poloidal component. Both sce- 
narios would happen on the Hall timescale 



THall 



4jm e eL 2 
cB 



(10) 



where n e is the electron density and e the elementary charge. For 
typical values in NSs we have T Ha ii ~ (10 4 -10 6 )^p yr. 

In situations of quasi-equilibrium, oscillatory modes with 
magnetic energy exchange between the different crustal field 



/ 2nt 

B (f) = B + 6B sin 

\ T"Hall 



and from Eq. ([8]) we have 



n = 3 - 4 — t c 
Bo 



3-8tt 



5B 1 2nt 
— cos 

^Hall Bo \T"Hall 



(ID 



(12) 



These oscillations are expected in young magnetars with 
Bo > 10 14 G. It causes corrections to the canonical braking 
index n = 3 of either positive or negative sign, according to 
Eq. (1121 1. These corrections are expected to be small because of 
their short r c , and to be increasingly important for objects with 
smaller P (old characteristic ages). In addition, as the star cools 
down, the drop in electrical resistivity may activate the Hall term 
even for normal pulsars Bo > 10 12 and the occurrence of oscilla- 
tory modes during a second Hall stage at late times is a natural 
outcome. First estimates suggest that 5B/Bq could be as large as 
~ 10 _1 . When r c » rnaii, the second term in Eq. (fT2l dominates 
and the magnetic field oscillatory modes should result in equally 
probable positive and negative braking indices with high abso- 
lute values. However, this needs to be confirmed by long term 
realistic simulations. 



4. Discussion. 

4.1. Magnetic field evolution in the pulsar population? 

From the whole popul ation in the ATNF Pulsar Catalogue 
(Manches ter et alj |2005), we extracted a sample of pulsars for 
which v has a quoted error smaller than 10%. From this pre- 
selected sample, we excluded all pulsars in binaries and those 
with very short periods (P < 15 ms), likely to be recycled, and 
we doubled check ed our list with the more recent review by 
Hobb s et alj d2010). the first large-scale analysis of pulsar timing 
noise over timescales > 10 yr, which led us to remove several 
more pulsars whose revised values were inconsistent with the 
ATNF data, or with larger errors. After this selection, our sam- 
ple contains 1 18 radio-pulsars, about half of which have negative 
braking indices. 

We show in the top panel of Fig.Q]the observed distribution 
of braking indices as a function of r c . We considered the char- 
acteristic magnetic field evolution timescale, Tg = 4r c /(3 - n) 
(see Eq. (O) and separated the sample into two groups, those 
with positive (middle pan el) and negative (bottom panel) tb- The 
eight pulsars discussed in Esp inoza et alJ d201 ll) are marked with 
squares. With all due cautions regarding the uncertainties asso- 
ciated to these measurements, it is worth mentioning some inter- 
esting trends visible in this plot: 

- All young objects seem to have always n < 3, which can be 
a hint of an increasing dipolar magnetic field. Q 

- All middle age objects (10 4 -10 5 yr), except Vela, that sur- 
vived to our conservative selection criteria have negative 
time derivatives of BqQ 



1 We use the values collected in Espin oza et al.1 ( 120111) for this plot. 

2 We did not include seven more objects, that fall in this region, with a 
quoted error of v smal ler than 10% in the ATNF pulsar database. They 
are not considered in Ho bbs et"aD i2010h . probably because they did 
not have 10 yr of accumulated data, although their values have not been 
reported to change. 
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Fig. 1. Top: braking indices (absolute value) as a function of t c 
for our sample of 1 18 pulsars. We represent pulsars with n < 3 
with crosses and pulsars with n > 3 with diamonds. Note that 
all objects marked with crosses and t c > 10 s yr have actu- 
ally n < 0. The 8 youngest objects discussed in Esp inoza et alj 
(1201 ll) are marked with squares. Middle: magnetic field evolu- 
tion timescale, tb = 4t c /(3 - «) for pulsars with n < 3. Bottom: 
Absolute value of (they are all negative) for pulsars with 
n > 3. 



- For old objects, there is no correlation at all, and there are 
similar numbers of objects with positive and negative deriva- 
tives of the field. The typical evolution timescales are in the 
range 10 3 -10 5 yr. 

The strong correlation between \n\ and r c seen in the old ob- 
jects of the top panel simply refiectis the definition of n, Eq. ([8]), 
with the additional piece of information that tb does not seem to 
be correlated with r c (see middle and bottom panels). 

Note that our criteria automatically select the objects with a 
dominant contribution of v in the ti ming phase re s iduals (i.e., 
the cleanest cubic lines in Fig. 3 of iHobbs et alJ {2010), e.g. 
B01 14+58). This introduces a bias toward objects with high 
braking indices, and against pulsars with residuals dominated by 
higher order terms (e.g. B0136+57) or quasi-periodic terms (e.g. 



B 164 2-03 or B 1826- 17) in the time-dependent phase dLvne et al.l 
I2010h . In the latter cases, the estimates of v are subject to larger 
uncert ainties, and n stron gly depends on the analyzed time in- 
terval dHobbsetalJ l2010). However, we have checked that the 
general trends do not change if the sample is enlarged by includ- 
ing pulsars with quoted errors in v up to 80%. This simply in- 
creases the statistics (> 300 sources) and includes some objects 
with lower value of \n - 3| (shorter \t b \). 



4.2. Expected evolution from theoretical models. 

We performed a series of numerical simulat ions with the 
magne t o-thermal evolutio n co de presented in | Aguilera et al.l 
d2008l) : iPons et alj d2009l) . and IVigano et all d2012l) , which in- 
cludes all relevant mi crophysical processes and follows the com- 
plex feedback between the physical mechanisms presented in a 
simplified way in the previous section. In addition to Ohmic dis- 
sipation and Hall drift effects, which are consistently included in 
the simulations, we assumed that all models suffer an episode of 
hypercritical accretion during their formation. 

We plot in Fig.|2]the evolutionary tracks of some theoretical 
models in the P-P diagram, compared to the pulsars of our sam- 
ple. We show results for three initial magnetic field configura- 
tions. Two of them are purely dipolar, although during the evo- 
lution other multipoles and toroidal field are naturally created. 
They differ by the magnetic field strengths: Bq = 10 13 G (model 
A, dashed lines) and 10 14 G (model B, dotted lines). In each case 
(models A and B) we compare results with three values of the to- 
tal accreted mass. The third initial configuration (model C, solid 
line) has an initial dipolar field of Bq — 10 13 G, as model A, but 
with an additional octupolar component of strength 3.5 x 10 14 G 
(at the pole). This model serves to make explicit the effect of the 
complex initial geometries. Assuming vacuum external bound- 
ary conditions, the toroidal field has to vanish at the star surface, 
but the dipolar component is still coupled to higher order multi- 
poles, which can also interchange energy with the dipolar mode. 

For the eight youngest pulsars, we indicate their predicted 
movement for th e next 2.5 kyr wi t h arro ws, assuming the present 
value of n, from lEspinoza et alJ (1201 ll) . remains constant. The 
direction of the tangent vector to a given track is related to the 
braking index. For reference, we indicate in the legend on the 
upper right corner that n — 3 and n — 1 imply constant inferred 
Bo, and constant characteristic age, respectively. 

The initial period P,„, which was assumed to be 0.01 s in 
all cases, only affects the early stage, while P < Pi„/t (with t 
being the real age), and P P,„. Models with deep submergence 
of the field but with different initial periods also have vertical 
trajectories shifted to the left/right depending on P,„, and quickly 
cross the range of P where the bulk of pulsars lie. At late times, 
tracks coming from the same model but with different P m are 
indistinguishable, typically converging after t ~ 30(q-^ b — ) 
yr, with B m being the initial magnetic field. Therefore, P,„ has 
an appreciable long-term effect only in the deep submergence 
case (M fl > 10~ 3 M o ), for which B m is strongly reduced: in these 
models, the first few 10 4 yr are spent in the vertical trajectories, 
with n <k 0. Note that in this reemergence phase there is no 
correlation between the real and characteristic ages. 

When reemergence of the field has almost been completed, 
the trajectories reach the high-P region (i.e., largest Bo) and 
progressively bend. The mor e extreme braking indi ces of PSR 
J1734-3333 (n = 0.9 + 0.2. lEspinoza et al.l fcoilh ) and PSR 
J0537-6910(n ~ -1.5. lMiddleditch et al.1 (120061) ) would be con- 
sistent with the last stage of the reemergence after a deep sub- 
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Fig. 2. Evolutionary tracks in the P-P diagram during the first 
3 Myr for: model A (dashed) and model B (dotted), both for 
different values of accreted mass M a = [0.1,5, 10] x 1O~ 4 M ; 
model C (solid) without accretion. 
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Fig. 3. Braking index (n) as a function of age for model A 
(dashed line), model B with M a = 1O~ 3 M (thin solid line), and 
model C (thick solid line). 



mergence into the inner crust. On the other hand, in the shal- 
low submergence models (M a ~ 1O -S -1O~ 4 M ), after accretion 
stage Bq =a Bi„. These tracks initially run almost along the iso- 
magnetic lines: pulsars with n slightly less than 3 are compatible 
with this scenario. 

Independently of the early reemergence phase (if any), tracks 
with the same B, n converge at middle-age, and have slopes cor- 
responding to n > 3, characteristic of the slow Ohmic dissipa- 
tion. During this phase, there is a correlation between real and 
characteristic ages, with typically t c being a factor of few longer 
than the real age. Some tracks show visible oscillations produced 



by the Hall activity when the star is cold enough (f > 10 5 yr); 
in particular, model C (solid line) clearly shows that any com- 
plex initial geometry may have a distinct signature on the timing 
properties of pulsars. 

We emphasize that we did not attempt to fit individual ob- 
jects: our purpose with this sample of models is simply to show 
that, with reasonable assumptions, it is possible to explain the 
variability in the observed range of braking indices of young pul- 
sars and to predict their evolutionary paths in the P-P diagram. 

4.3. Braking index and evolution timescale for realistic 
magnetic field evolution models. 

In Fig.[3]we plot the braking index evolution for three represen- 
tative models, Note that the horizontal axis in Fig. [TJ represents 
the characteristic age r c , while in Fig. [3] we show our results as 
a function of the real age of each model, so that a direct com- 
parison is not possible. We also used a linear scale in this plot to 
show the quasi-periodic oscillations during the long-term evolu- 
tion more clearly. For fields < 10 13 G and simple dipolar geome- 
tries (dashed line) the braking index at late times is n > 3, but 
its absolute value is low. In contrast, for strong dipolar fields, or 
for weak dipolar components but with strong higher order multi- 
poles, Hall-drift induced oscillations appear sooner or later and, 
in some situations, have large amplitudes that result in very high 
absolute values of the braking index. It is particularly interesting 
to compare models A and C (dashed and solid line), which have 
the same initial dipolar component, representative of a typical 
pulsar (at 10 6 yr the dipolar field is about 3-5 x 10 12 G). However, 
the presence of a strong octupolar component at birth results in a 
radically different braking index behaviour, even if P and P are 
similar. The amplitude of the oscillations and whether the modes 
are damped or excited depends on particular details of the small- 
scale structure of the magnetic field, which is unknown. The 
exact age at which these oscillatory modes are excited is con- 
nected to the temperature of the star, and therefore to its cooling 
history and internal physics (neutrino emission processes, super- 
fluid gaps, etc.). For the standard cooling scenario this happens 
at w 10 5 yr. We can also observe that as the star evolves and the 
magnetic diffusivity decreases, the frequency of the dominant 
mode may vary (see solid line). For example, at early times in 
model C we observe variability on shorter timescales (few kyr), 
while later we see longer oscillation periods (10 5 yr). After 1.5 
Myr higher frequency modes seem to be growing again. 

4.4. Timing residuals. 

In timing analysis of radio pulsars, an important piece of infor- 
mation is the study of the residuals, i.e, the phase differences 
between the observed signal and the best-fit model, including 
the frequency and frequency derivative. As shown in detail in 
iHobbs et al J (12010). there is a rich variety of shapes in the resid- 
uals of radiopulsars. As we have shown, the strong Hall-induced 
interplay between different multipoles and the toroidal field pro- 
duces a complex evolution of Bo- To compare with observational 
timing analysis, we proceeded as follows: first, we choose a short 
interval of 40 yr in our simulations, centered on a fixed time de- 
noted by fo- This interval is similar to the longest periods for 
which phase-coherent timing analysis for radio pulsars can be 
performed. Then, from our theoretical Bo(t), we obtain v(t) by 
integration of the classical spindown formula (Eq. [6]) assuming 
an orthogonal rotator, and the phase <t>(?) = f v(t')df. Finally, 

" Jin 
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Fig. 4. Phase residuals for model C around fo = 10.93 kyr (upper 
panel) and to = 100.5 kyr (lower panel), with the corresponding 
value of n indicated in each case. We show residuals after re- 
moving the best-fit solution including quadratic terms (solid) or 
including cubic terms (long dashes); dotted lines show the cubic 
residuals obtained assuming a constant Bo- 

we fit our synthetic time-dependent phase with a quadratic func- 
tion 

1 9 

O(0 = O + v f (t - t ) + -v f (t - to) 2 . (13) 

The results of the fit, v/ and Vf, are of course in the range of 
values of the real quantities in the time interval. 

In Fig. |4] we show the phase residuals for model C (solid 
lines) at to = 10.93 kyr (top) and f = 100.5 kyr (bottom). During 
the two analyzed intervals, the mean values of the dipolar mag- 
netic field are B = 8.1 xlO 12 G and B = 5.3xl0 12 G. The cubic 
shape indicates that the residuals are dominated by the next term 
in the Taylor series (i.e., red noise). Including in the fitting func- 
tion a cubic term allows one to measure v and the braking index 
(n = 0.39 and n = 3.45, respectively). With dashes we show 
the fourth-order residuals after subtracting the third-order term 
in the fitting function. For comparison, we also show with dot- 
ted lines the third-order residuals obtained assuming a constant 
value Bo(t) = Bq (that leads to n = 3). 

In the first time interval, the change of magnetic field, SBo = 
1.3 x 10 10 G = 1.6 x 10~ 3 B(), is strong enough for the cubic 
residuals to be visibly different from the constant field case. The 
value of |v| (and of \n\) is low, and the residuals of fourth order 
are significant. In the second time interval, 6Bo = 1.4 x 10 8 G = 
2.6 x 10~ 5 Z?o, and n is close to 3. As a consequence, the deviation 
from the constant field spin-down behavior is slow and fourth- 
order residuals are orders of magnitude weaker (dashed line). 

To investigate the effect of short-time irregularities, we re- 
peated the process but artificially added two additional sources 




-20 -10 10 20 
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Fig. 5. Comparison of the cubic residuals of the theoretical 
model (solid line) with the result after adding a sinusoidal per- 
turbation with frequency 0.1 yr (dashed) or a random pertur- 
bation (dot-dashed). For both perturbation forms we considered 
a maximum amplitude of(8B /2) = 6.5 x 10 9 G. 

of 'noise" (Bp) to our theoretical values of Bo(t), with amplitude 
of the same order as the theoretical variations 6Bo- In Fig. [5] we 
show the residuals obtained by adding a sinusoidal perturbation, 
B p = (6Bo/2) sin(27rf/r), with T — 10 yr (dashed lines), or a 
random perturbation of maximum absolute value 6Bq/2 (dash- 
dotted lines). In our analysis, we sampled values every ~ 10~ 2 
yr, comparable with the typical integration time in observational 
data. We checked that increasing the time interval strongly re- 
duces the contribution of the random noise at a fixed amplitude, 
as expected. Similarly, the periodic short-term noise is reduced 
when the integration interval becomes much longer than the per- 
turbation period. 

Note that small perturbations of Bo produce visible effects, 
but the cubic residual still dominates. In constrast, in many pul- 
sars a smooth cubic phase residual is clearly seen, which means 
that the spin-down is very stable, short-term irregularities are 
completely negligible (or averaged out in the long observation 
period), and the measure of a braking index (even a high value) 
dominated by the secular evolution is robust. Glitches should 
ap pear as cusp s in the se plots (e.g. see data for PSR B0 154+61 
in lHobbs et alJ (120101) '). but a more detailed investigation in this 
direction is beyond the scope of this paper. 

5. Conclusions. 

According to our current understanding of the magneto-thermal 
evolution of neutron stars, the time variation of the magnetic 
field notably affects the braking index of pulsars and contributes 
to the so-called red noise in the timing residuals. The qualitative 
picture seems to agree with the data: i) there is a short stage in 
which the dipolar field appears to be increasing, possibly reflect- 
ing the reemergence of magnetic field after initial accretion; ii) 
a stage dominated by Ohmic dissipation, while the star is still 
warm, and during which we always expect n > 3; iii) once the 
star has cooled down, after about 10 5 yr, oscillatory Hall-drift 
modes in the crust are excited, and quasi-periodic oscillations of 
the braking index are consequently expected. 

This is by no means the only contribution to the timing noise, 
nor necessarily the dominant one. For example, changes in the 
effective moment of inertia as the superconducting region of the 
core grows and unknow n mechanisms for alignment of the rota- 
tion a nd magnetic field (Tauris & Manchester 1998; Young ^tai] 
2010) would have similar long-term effects. In addition to this 
secular evolution, short-term magnetospheric effects, such as os- 
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dilations of the inclination angle around the equilibrium po- 

i 1 1 1 — L -ii — i 

sition, or magnetospheric changes (ILi et al. 2012a b), could be 
also responsible for the strong q uasi-periodic fe atures seen in 
phase residuals for several cases dLvne et all 2010). 

Our results show that realistic magnetic field evolution mod- 
els predict temporal variations consistent with the observed 
trends. In particular, configurations similar to our model C could 
explain the wide range in positive and negative values that are 
observed. This model has a moderate dipolar component (10 13 G 
at birth, which becomes a few times 10 12 G at middle age), typ- 
ical of pulsars, but a strong octupolar component, large enough 
to activate Hall modes. The same generic behavior is expected 
for other models with strong, smaller scale components, even if 
the dipole is weak. However, it must be mentioned that the mod- 
els we studied predict high values of the braking index (10-100) 
but not extremely high values (10 3 -10 4 ), as reported for pulsars 
with old characteristic ages. The absolute value of the braking 
index at late times, which can be of either sign, is not simple to 
predict, since it depends on particular details of important com- 
ponents of the magnetic field (higher order multipoles, toroidal 
component). Perhaps different field geometries (e.g., without ax- 
ial symmetry) and/or magnetospheric corrections may account 
for the discrepancy. 

We emphasize that there is absolutely no reason to expect 
that the magnetic field remains constant during a pulsar lifetime, 
and we must abandon the oversimplified models that assume 
constant (in time) pure dipolar fields. On the other hand, since 
there is a strong interplay between the magnetic, thermal, and 
rotational evolution of neutron stars, there are unexpected ways 
to obtain information. If future observations (or a deeper analy- 
sis of archival data) increase the number of reliable measures 
of braking indices, and we can firmly establish at which age 
the transition from the purely diffusive to the oscillatory regime 
happens, there is coded information about the cooling history of 
neutron stars of potential interest to place constraints on the in- 
terior physics. In this respect, it is also important to mention that 
precise phase-connected timing solutions can also be performed 
in other bands, as has b een recently done for 17 y-ray-selected 
pulsars dRav et al.ll201 lb . using the Large Area Telescope (LAT) 
on the Fermi Gamma-ray Space Telescope. This and other future 
high-energy missions can enlarge our data sample and improve 
the precision of our measurements. 
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